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We analyse the thermodynamic quantities of some simple oscillator systems like bosonic harmonic 
oscillator, fermionic harmonic oscillator and a supersymmetric harmonic oscillator (which is a combi- 
nation of a bosonic and a fermionic oscillators), and discuss in detail about the nature of the specific 
heat, internal energy and entropy of these systems both at low and high temperatures. Also we 
have studied the behavior of the thermodynamic properties of the well known Landau-Fock-Darwin 
(LFD) problem viz., an electron in the combined presence of a uniform magnetic field acting perpen- 
dicular to the plane of motion of the electron and an isotropic and cylindrically symmetric parabolic 
potential well in the directions normal to the field with a natural confinement frequency wo- We 
have succinctly derived the respective thermodynamic quantities of LFD model with and without 
the intrinsic spin angular momentum of the electron. For the sake of simplicity and convenience we 
have taken the simple case of an electron carrying spin ~. 
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I. INTRODUCTION 

Thermodynamics; the phenomenological branch of 
physics, is a summary of the properties that real phys- 
ical systems (in thermal equilibrium) exhibit. Thermo- 
dynamics was established in the 19th century by a group 
of pioneers with the formulation of the law of conserva- 
tion of energy relating the equivalence of work and heat 
as stated by Robert Mayer and Hermann HclmholtzpQ. 
The total dynamics of a system can be defined indepen- 
dantly in terms of the main constituents of thermody- 
namics like entropy, temperature and the three laws con- 
necting these state variables. This could be a reasonable 
argument for the robustness of thermodynamics over the 
last century, eventhough some of the profound peers like 
classical mechanics and electrodynamics underwent a lot 
of structural changes. 

Thermodynamics shows its full glory when it couples 
with Statistical physics, particularly for a microsopic sys- 
tem. Statistical physics is that part of physics which de- 
rives emergent properties of macroscopic matter from the 
atomic structure and the microscopic dynamics. Emer- 
gent properties here means those properties like temper- 
ature, perssure, dielectric and magnetic constants etc, 
which are essentially determined by the interaction of 
many particles (atoms or molecules). Exclusively, these 
emergent propeties are typical for many-body systems 
and they do not exist (in general) for microscopic sys- 
tems. Strictly speaking, statistical physics is the bridge 
between the microscopic and the macroscopic world and 
it provides methods for calculating the macroscopic prop- 
erties (for eg., like the specific heat), from the micro- 
scopic information( like the interaction energy between 
the particles) . Statistical mechanics introduces probabil- 
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ities into physics and connects them with the fundamen- 
tal physical quantity 'entropy' (The name 'entropy' was 
given after the work of Clausius in 1865 [5] and also James 
Clerk Maxwell [3] in 1878 through his " Gedankenexper- 
iment" . In fact, Maxwell proposed the name "Statisti- 
cal Mechanics".). The coalescence of thermodynamics 
and statistical physics led to the emergence of Statisti- 
cal thermodynamics, established through the well known 
Boltzmann-Planck [4] formula for the Clausius entropy 
S = kg In W (W is the thermodynamic probability) and 
the ergodic hypothesis by Boltzmann [S], also indepen- 
dantly by Josiah Williard Gibbs [6] , the father of the en- 
semble theory. In classical statistical mechanics we did 
not know the exact microscopic state of the system (lo- 
cation in the phase space) . We made use of the ergodic 
hypothesis and replaced the time average of a physical 
quantity with an ensemble average ie., an average over 
many equivalent systems. 

All these developments in classical statistical mechan- 
ics almost coincided with the revolutionary concept of 
quantum theory by Max Planck, that changed the entire 
structure of physics. It would have been surprising if sta- 
tistical mechanics could have escaped the repercussions 
of the quantum revolution. But the whole structure of 
statistical mechanics was overhauled by the introduction 
of the concept of indistinguishability of (identical) par- 
ticles. But the concept of ergodicity remains untouched 
in the quantum case also. Microcanonical and canonical 
ensembles retained their positions with some modifica- 
tions and grand canonical ensemble made a debut. The 
statistical aspect which already exists in classical statis- 
tical mechanics in view of the large number of particles 
present in the system, has been augmented by the sta- 
tistical aspect came from the probabilistic nature of the 
wave mechanical description. All these new concepts re- 
formulated the ensemble theory with the introduction of 
density matrix, which is the quantum mechanical ana- 
logue of the density function of the classical phase space, 
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first introduced by Landau [7] and later by von Neumann 
[8], later on more rigorously by Dirac [9]. 

With these remarks we are led to analyze the quantum 
thermodynamic properties of some simple oscillator sys- 
tems like a bosonic harmonic oscillator, a fermionic har- 
monic oscillator and finally a supersymmetric harmonic 
oscillator (which combines a bosonic and a fermionic os- 
cillator). Harmonic oscillator is the model system of 
model systems and it is effectively a single particle prob- 
lem, ie., it is a problem where the particle is bounded 
within a quadratic potential. Hence it is always intutive 
to study the basic quantum thermodynamical properties 
of a particle trapped in a quadratic potential irrespec- 
tive of the inherent character of the particle (Boson or 
Fermion). Classically speaking, the harmonic oscillator 
system is like a particle moving in the x-direction con- 
nected by a spring to a fixed point. Its potential en- 
ergy is V(x) — ^|-, where k is the force constant of the 
spring attached to the particle. In the case of a clas- 
sical oscillator, the amplitude A of the periodic motion 
can take a continuum of values, and hence the energy 
E = ^muj 2 A 2 , starting from E = 0. But quantum me- 
chanical harmonic oscillator can only take certain dis- 
crete energy values since the energy eigen spectrum is 
fully quantized or discretized as we will see in the forth- 
coming sections in detail. 

Our main aim is to study the statistical thermodynam- 
ical quantities like internal energy, entropy and specific 
heat (an attribute of the third law of thermodynamics 
|10j). Third law says that the specific heat and the lim- 
iting entropy of a quantum system vanishes as the equi- 
librium temperature is approa- 
ched. Perhaps Einstein's work in 1907 [TT] (considered as 
the origin of quantum statistics) proposed that quantum 
effects lead to the vanishing of the specific heat at zero 
temperature. In addition to that, we consider the effect 
of magnetic field on the thermodynamic properties of an 
electron in the combined presence of a uniform magnetic 
field and a parabolic potential. 

The quantum mechanical problem of an electron in a 
magnetic field was first solved by Landau[H]. Peierls[T3"] 
elucidated the concept of diamagnctism of free elec- 
trons under the influence of a magnetic field, followed 
by Darwin [i"4j who explained the role of boundary for 
the recovery of correct diamagnetism. Even from the 
classical mechanical point of view this problem was quite 
interesting and surprising. The famous van Leeuwen|15j 
theorem was a fatal blow to some of the older theo- 
ries of diamagnetism that abounds with subtle pitfals 
because of the reason that the correct bulk susceptibil- 
ity calculation needed a boundary effect. Later on the 
problem of an electron in a magnetic field has been im- 
mensively exploited and has got innumerable attention 
in terms of quantum transport phenomena like quantum 
hall effect [16]. 

Another paradigm, an electron in the combined pres- 
ence of a uniform magnetic field and a parabolic potential 
was comprehensively solved by Fock|17) (nowonwards we 



call this model as Landau- Fock-Darwin (LFD) model). 
More recently LFD model has been extensively studied 
under the influence of a dissipative quantum heat bath 
of non-interacting oscillators for the calculation of dissi- 
pative Landau diamagnetism [15] and to verify the third 
law of thermodynamics [19]. In this paper we are not in- 
terested in the much studied Landau diamagnetism, but 
we will be looking at the equilibrium thermodynamics of 
the LFD problem. 

We restrict ourselves to equilibrium statistical mechan- 
ics, which aims at characterising the equilibrium state in 
terms of given constraints, without regard to the mech- 
anism relating to how this equilibrium state is brought 
about. Since the systems we are interested in are those 
with the number of particles fixed, we can use canonical 
ensemble formalism to obtain the thermodynamic quan- 
tities. 

With the preceding introduction the organization of 
the paper is as follows: In Sec. 2, we review the well known 
thermodynamic quantities which we want to look at. In 
the following subsection, we revisit the quantum thermo- 
dynamics of a simple bosonic harmonic oscillator. Other 
two subsections mainly decribe the thermodynamics of a 
fermionic harmonic oscillator and a supersymmetric har- 
monic oscillator. In Sec. 3, we analyse the problem of 
an electron in the combined presence of a magnetic field 
and a parabolic confinement potential. Here we discuss 
the two cases viz., electron with and without its inherent 
spin. In Sec. 4, we discuss the results and finally, Sec. 5 is 
devoted to the concluding remarks. 



II. STATISTICAL THERMODYNAMIC 
QUANTITIES REVISITED 

Canonical ensemble description is the most useful one, 
as in most practical cases it is possible to control the 
temperature and not the energy of the system. Canoni- 
cal ensemble theory depicts the situation where one has a 
system with a fixed number of particles, but can exchange 
energy with the environment with which it is attached. 
More precisely, a system is kept at a fixed temperature 
when it is in thermal equilibrium with a much larger sys- 
tem, referred to as a heat bath. Clearly, the energy of the 
system under study is not a constant, but it has a well 
defined average energy and the magnitude of the heat 
energy exchanged with the bath is always small com- 
pared to the average energy. We thus have to determine 
probabilities for states with different energies. The rel- 
ative recurrence rate of states with different energies is 
given by their Boltzmann factor, namely, by the quantity 
e -E/k B T ^ w ]j ere 25 is the energy of the small system, and 
ks is the Boltzmann constant. In the canonical ensemble 
the system is assigned a temperature, and any state can 
appear in it. With the small description, we define an all 
important concept in the canonical ensemble, an ensem- 
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ble averaged quantity known as the partition function 



Z = 



-^^(microscopic state) 



all microscopic states 



(1) 



Partition function as such does not have a physical 
meaning, but this eventually clarifies the connection with 
thermodynamics. Once we obtain the partition function, 
then it is easy to derive the Free energy which is given 
by 



F = -k B T\nZ 



(2) 



From this formula we can derive all the thermodynamic 
properties, also the Internal energy (or average energy) 
can be written as 



U 

The Entropy 

S = k B {lnZ 
and finally Specific heat 



dp 



InZ 



d/3 



InZ} 



,dU 
~d~p ' 



Or 



8S_ 

dp 



(3) 



(4) 



(5) 



(6) 



Now we have the required formulae for the calculation of 
the thermodynamic quantities of the systems which we 
are interested in. 



A. Simple (Bosonic) Harmonic oscillator 

Let us briefly describe the well known quantum har- 
monic oscillator. We have a particle of mass m, trapped 
in a parabolic well and the entire arrangement vibrates 
with a natural frequency u. The whole system is in equi- 
librium with a thermal reservoir at temperature T. The 
harmonic oscillator is a particular case in which the de- 
gree of freedom takes discrete values, but the number of 
different values are infinite. The states of the harmonic 
oscillator are all non-degenerate. Probability of a state 
with energy E n is same as the probability of the oscilla- 
tor having energy E n , with n as the oscillator quantum 
number. The probability can be defined as 



P(E n 



1 



exp(-/3£' n ) 



where the partition function can be calculated as 
Z = Tr{e- mB ) , 



(7) 



(8) 



Where % B denotes the Hamiltonian for the bosonic har- 
monic oscillator (The subscript B is for bosonic). We 
can write down the Hamiltonian of the bosonic harmonic 
oscillator in the (x,p) space and is given by 



P 

2m 



1 



2 2 
m id X 



(9) 



It is also possible to write down the Hamiltonian in 
terms of the creation and annihilation operators a B and 
a b (Mathematical details can be found in some standard 
texbooks like Quantum Mechanics by Zettili [5D]) which 
is given by 

U B = fw(a B a B + -) = hu{N B + ~) , (10) 



The subscript B is to identify the operators as bosonic 
We know a B \n B ) = \fn~B 
y/nE\n B - 1 
as 



l\n B + 1) and a B \n B ) = 
where the wave function can be written 



n B 



7=T( a B)™ s l >s- Also a B a B = N B , where N B 
is the number operator. The Bosonic operators obey 
the commutation relation [as,a^]_ = 1. We know that 
Ob|0) = 0. The Hamiltonian operating on the state \n B ) 
can be written as H B \n B ) = (n B + ^)hui\n B ) . It is ob- 
served that the energy-eigen values are E n = (n B + ^)huj. 
Note here that the energy eigenvalues E n are positive, 
just as for a classical harmonic oscillator. But in the 
quantum case, the energies of the oscillator are quan- 
tized in intervals of hu>, starting from a non-zero value 
E = ^hu>, the zero point energy. 

With all these requirements, the partition function can 
be calculated as (cf., Eq.(8)) 



Z = Tr(e- mB ) 

riB—oo 

= cxp(— [3huj(n B 

n B =0 

1 



2sinh(^) 



(11) 



From the partition function we can derive all the ther- 
modynamic quantities, and are given by: 
Free energy 



F = 
Internal energy 

T 4 



~2 



k B T\n{l~e- phw ) , 



U 



Entropy 



S 



(ihuj 



hu phuj 

— cothl 

2 v 2 



k B eP riuJ - 1 
alternatively, 



-ln(l 



), 



S Phw ,,,PT>Ws , ro . ,,P^s., 
- = — coth( — )-ln[2smh( — )] 



(12) 



(13) 



(14) 



(15) 
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and finally the specific heat 



Now the Hamiltonian can be written as 



Cv 
k B 



(/^) 2 cosech 2 (^) 



a f3huj 



(16) 



Thus we have derived all the thermodynamic attributes 
of a simple bosonic harmonic oscillator. 



H F = hw(a F a F - ~) = hw(N F - ~) . (21) 

Now if we assume \n F ) as an eigenstate of N F , then 
we have Nf\hf) = n F \n F ), np — 0, 1. Let the ground 
state with no quantum be denoted by |0) and it satisfies 



7V F |0) = then =► H F \0) = 



(22) 



B. Fermionic Harmonic Oscillator 

We have described the bosonic harmonic oscillator in 
one dimension with a natural frequency w. The Hamilto- 
nian of such an oscillator can be written in terms of the 
creation (a B ) and annihilation (ag) operator that obey 
a commutation relation. The Hamiltonian has a sym- 
metric form since we are dealing with Bose particles, so 
the states must also have a symmetric form. The opera- 
tors corresponding to bosons obey commutaion relation, 
and that for fermions satisfy anti-commutation relations. 
That means the fermionic system have an inherent anti- 
symmetry associated with it. Thcrfore let us try out a 
Hamiltonian for the fermionic harmonic oscillator with a 
frequency uj of the form |21j 



u ^ i t 
Hf = -^{a'pap 



a F 



4) 



(17) 



The subscript F is for fermionic. Now, it is interesting 
to see that, the fermionic Hamiltonian Tip = —^f> if we 
assume the operators ap and a F obey the same commu- 
tation relations [af ,a^]- = 1 as bosons. Then there is 
no dynamics associated with it. So let us assume that 
the operators satisfy anti-commutation relations, 



[op, op] 



[a F ,a F ] + = 1 = [a F ,a F ] + 



[a F , a F } + , 



(18) 



Anti commutators are by definition symmetric. It is 
quite obvious from the anti-commutation relations that, 
in such systems the particles must obey Fermi - Dirac 
statistics. The justification is as follows. If we identify 
a p and a F as annihilation and creation operators for such 
a system, then immediately we can define a number op- 
erator 



a/pCLF 



(19) 



then it is easy to see, Np = Np or Np(Np — 1) = 0. 
Therefore, the eigenvalues of the number operator can 
only be zero or one. This is the reflection of the Pauli 
principle or the Fermi-Dirac statistics. Also it proves that 
the anti-commutation relations are the natural choice for 
a fermionic system. We also have 



[of, N f } + = a F 



p ■ Np\ + 



(20) 



simiarly, the state with one quantum is identified to be 
|1) and it satisfies 



N F \1) = |1) then 



Up\l) 



ID 



(23) 



we have a F \0) =0 



|0) 



ID 



|1) =0. Hence, in 



this case the Hilbert space is two dimensional. Also, one 
more interesting point is that, for a fermionic harmonic 
oscillator, the ground state energy has the opposite sign 
from that of a bosonic harmonic oscillator. It is possible 
to write the Hamiltonian of a bosonic harmonic oscillator 
in the position and momentum coordinates, which basi- 
cally equals its classical counterpart. But fermions have 
no classical analogue, and as a result we cannot directly 
write down a Lagrangian for the fermionic oscillator with 
the notions of coordinates and momenta. In order to do 
that, we need the notion of classical anti-commuting vari- 
ables. But such variables have been extensively studied 
in Mathematics under the name of Grasman variables 
|22j . We are not interested in explaining the details of 
the Grasman variables here. 

As we have identified, the Hamiltonian of a fermionic 
oscillator given by Eq.(21). The fermionic oscillator is 
like a two level system with energy eigenvalues 



~~2 



~2 



E, - ~ • E. ^ . (24) 
Now the partiton function is trivial and is given by 



Tre 



-mi 



2 cosh( 



(25) 



So forth the thermodynamic quantities: 

Free energy 

F = -i(ln2cosh(^)) 
Internal energy (or average energy) 



(26) 



U = — 



2 

hid 



2 tanh (^") 



Entropy 
S 



ln2cosh(— — ) 
k B 2 



(^— )tanh( 



(27) 



(28) 
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alternatively, 

-=ln( e ^ + l)-/3^-^— , (29) 
and finally the specific heat 

% = = (^) 2 p£^ ■ (30) 

The difference in specific heat expressions of the bosonic 
and the fermionic oscillator is just a matter of sign in the 
denominator. 

C. Supersymmetric Harmonic Oscillator 

Supersymmetric oscillator is a simple toy model in 
quantum field theory, and it is a combination of a bosonic 
and a fermionic oscillator with the same natural fre- 
quency oj 23 . Or in other words, it is a system where 
one boson and one fermion get trapped in a harmonic well 
and the whole arrangement vibrates with a frequency uj. 
The energy eigenstates correspond to a boson with indef- 
inite energy eigenstates and a fermion with two (either 
zero or one) energy eigenstates. The Hamiltonian is given 
by 

Hs = Hb + Hf = hu](aJ B aB + a F ap) • (31) 

with Nb — a B a B and Np — a F ap. Now Eq.(31) can be 
expressed as 

Hs = Huj{N b + N F ) . (32) 

It is clear from the above expression that the energy eigen 
states of the system will be the eigenstate of the number 
operators Nb and Np. Consequently, let us define 

\N B ,N F ) = \N B )® \N F ) , (33) 

where 

Nb\ub) = TIB | tib) , Ub = 0,1,2, 

N F \np) = np\np) , uf = 0, 1 . (34) 

Hence we may write 

Hs\n B ,np) = E nBtnF \n B ,np) 

= hhj(nB + np)\nB,np) . (35) 

So the supersymmetric harmonic oscillator has energy 
levels E n = hwn, n = ub + np. This resembles an 
ordinary harmonic oscillator, but without ground state 
energy. It is very much clear that the ground state energy 
of the supersymmetric oscillator vanishes, ie., 

£o,o = 0, (o B |0> = a F \0) = 0) . (36) 

In the realm of supersymmetry theories, vanishing of 
the ground state is not an unusual phenomena, it is a 



consequence of the supersymmetry of the system. From 
Eq.(35), except for the ground state, all other energy 
eigenstates are doubly degenerate. Precisely, the states 
Ins, 1) and |tib + 1,0) have the same energy for any 
value of ns- Or, there is one state with n F = and 
jib = n and another one with np = 1 and tib = 71 — 1, 
the same as an ordinary harmonic oscillator with ground 
state energy subtracted. So except for the vaccum state 
all states come in pairs of one fermionic and one bosonic 
state with same energy. Admitting to the supersymme- 
try theory, this degeneracy is also a consequence of the 
supersymmetry of the system. At this point we will dis- 
cuss a little bit about the supersymmetry of the system. 
In the theory, we have two more operators (fermionic) 
to describe (generally these operators are known as the 
supercharges) : 

Qi = a B a F and Q 2 = apa,B ■ (37) 

Since the bosonic operators commute with the fermions 

[Qi,H s \- = and [Q 2 ,H S ]- - . (38) 

Therfore Qi and Q2 define conserved quantities of the 
system and would correspond to the generators of sym- 
metries in the theory. Supercharges act on the paired 
states by transforming a bosonic state to a fermionic state 
and vice versa. The only exception is the ground state 
which is annihilated by the supercharges. This is a typ- 
ical situation for the supersymmetric system. There ex- 
ists models in which one has a supersymmetry algebra, 
but ground state has nonzero energy. In the language 
of field theory, the ground state is not supersymmetric, 
therfore the model is not supersymmetric, supersymme- 
try is spontaneously broken. The anti-commutator of the 
supercharges is given by 

[Qi,Q2}+ = ^Hs ■ (39) 

Certainly, Eqs (37), (38) and (39) indicates that, Qi, 
Q2 and Hs constitute an algebra which involves both 
commutators and anti commutators. This algebra is 
known as the graded Lie algebra. This algebra have a 
beautiful consequence, namely, if the ground state is in- 
variant under the supersymmetry transformations (ie., 
Oi|0) =0 = Q 2 |0) and 

(0\H S \0) = M0|QiQ 2 + Q2Q1IO) = 0), then the ground 
state vanishes in supersymmetry theory. Q 2 is really a 
hermitian conjugate of Q\. 

Energy eigenvalues of the supersymmetric oscillator is 
E n = hum, n — ns + np . (40) 
It is now straight forward to calculate the canonical par- 
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tition function: 



* = E 



-/3ftn 



n=0 



71=1 

e^ ri " + 1 



-phn 



(41) 



Or 



Z = coth(^). (42) 

Now as usual the thermodynamic quantities follows: 
Free energy 



1 , , .Bhu. 

F = ncoth(- — ) , 

Internal energy (or average energy) 

U hx 



smh( f3hu>) ' 



Entropy 
S 



,f3hw, 



= lncothf ) + (/3Hu})cosech(f3Huj) , 

k B 2 

and finally the specific heat 

C v 



k. 



(0hu>) cosech(j3hw) coth.(j3hu>) 



(43) 



(44) 



(45) 



(46) 



In an alternative form, specific heat can be written as 



CV , 2 e^ ( e 2 ^ + l x 



(47) 



This expression can be represented in terms of the 
bosonic specific heat as 



(%s = 2(% B (e^ + l) 
Kb k B 



(48) 



The specific heat has more bosonic character than 
fermionic. Hence the variation of the specific heat with 
temperature is expected to be more bosonic in nature. 



III. HARMONICALLY BOUND ELECTRON IN 
A MAGNETIC FIELD 

The problem of an electron in a uniform and homo- 
geneous magnetic field attained tremendous attention 
nowadays in mesoscopic physics [M| and quantum dots 
or the quantum confined nanostructures [25] . It was Lan- 
dau who emulated this problem and conveyed that the 
problem of an elctron in a magnetic field is essentially 
quantum mechanical and can be solved exactly through 



the first principles of quantum mechanics. We are in- 
terested in the thermodynamics of this problem with an 
additional harmonic oscillator potential. 

Since the energy eigen spectrum of the free electron un- 
der the magnetic field (both in the absence and presence 
of the parabolic confinement) shows simple harmonic os- 
cillator structure, we must commit to say that the quan- 
tum mechanical formulation is identical to that of a har- 
monic oscillator. Moreover, harmonic oscillator potential 
is the only case where simple results can be obtained. It 
allows one to discuss in detail the effects of the boundary 
and the way in which the bulk limit is attained. This 
is the reason we are interested in deriving the thermo- 
dynamic properties of this Landau-Fock-Darwin model 
(LFD). It is worth to study the effect of spin of the elec- 
tron trapped, on the model, so we analyse two cases viz., 
the LFD model with and without spin. For the sake of 
simplicity we have considered the case where the electron 
possesses spin ~. 



A. Without Spin 

Here, the magnetic behavior of the electron has been 
considered. Permanent magnetic dipole moment of elec- 
tron is ignored and only the influence of the external 
magnetic field on the electron is taken into account. The 
Hamiltonian for a non-relativistic electron of mass m and 
charge e confined in a two dimensional isotropic harmonic 
bowl (parabolic well) of frequency loq can be written as 



rLr, 



2m 



(P 



6 Af 

c 



22 
-mw„r 



(49) 



Where P and r are two-dimensional vectors and A is the 
magnetic vector potential. Using the "symmetric gauge" 
A = (Hy, —Hx,0), we can write the Hamiltonian in 
Eq.(49) as 

Ho 



1 



[(Px 



2c ' 



2m' 

■ 1 2 / 2 , 2 \ 

+ o^qC 35 +y > ■ 



(Pv 



exH s 

~2c~' 



(50) 



The energy eigen values of this particular problem is 
known as the Fock-Darwin spectrum which basically has 
a form 



= hJuo + (y) 2 (™i +n 2 + l) + ^hw c (ni - n 2 ) . 

(51) 

Where oj c = — , the cyclotron frequency, c is the veloc- 
ity of light. The Fock-Darwin energy spectrum is closely 
related to the energy eigen spectrum of an ideal rotat- 
ing BEC, with the trap is circular in cross section. This 
spectrum is interesting because it shows how smoothly 
the spectrum interpolates between the bound states of 
a two-dimensional harmonic oscillator (when lo c — 0) to 
the case of the Landau levels of a free particle in a mag- 
netic field (when ujq — 0). The Hamiltonian (cf., Eq. 
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(50)) can be effectively canonically transformed into the 
sum of two independent one dimensinal harmonic oscil- 
lators with two frequencies uj + and w_ either by intro- 
ducing the generating function technique of Valatin[2l)] 
or by quantum mechanical operators [27] • We define two 
operators 



1 



(52) 



with I = \ — is the Fock-Darwin length. Now it is easy 
to write the Hamiltonian in terms of these operators and 
it turns out to be 

Wo = (a i + a + + -)huj + + {a)_a_ + ~)tkV- , (53) 



where 



2 ' 

Y ' 



(54) 



finally Specific heat 



Cv 
k B 



cosech 2 ( ) 



cosech 2 ( 



2 



) 



(60) 



In the limit of the vanishing magnetic field (cj c = 0), 
this problem reduces to the problem of two isotropic har- 
monic oscillator of same frequency. Also in the limit of 
vanishing oscillator potential (wo — > 0), we can recover 
the the thermodynamics of the original Landau problem 
(ie., the electron in a magnetic field). In the vanishing 
oscillator potential limit (ie., loq —> or w_ — > 0) we can- 
not recover the partition function of the original Landau 
problem form Eq.(56) because Eq.(56) diverges in this 
limit, since the system becomes translationally invari- 
ant in space. Recovery of partition function for Landau 
problem is well expalined in 19J. Apart from the parti- 
tion function none of the thermodynamic quantities are 
plagued by this issue. 



with u> 
,t 



U ^ L . Now identifying ri\ — a\a + and 



ri2 = a_ci- gives rise to the energy eigen values 



E ni ,n 2 = (n-i + + (n 2 + 7^)^- 



(55) 



which is the same as Eq.(51). In the absence of a 
magnetic field the Fock-Darwin levels are degenerate, 
uj + = u_ = luq, whereas a strong magnetic field leads 
to the formation of the structure of Landau energy lev- 
els, represented by the cyclotron energy hw+ « Huj c . 

Thermodynamics of this problem has been studied by 
Jishad et al [19], where they have used the path integral 
approach to calculate the partition function. We obtain 
the partition function in the canonical ensemble as 



Z = Z+Z- = 



1 



1 



2s inh(^±)'2sinh(^) 



Now the thermodynamic quantities: Free energy 



(56) 



B. With Spin 

Assume that the electron possesses an intrinsic spin of 
value \ha and a magnetic moment [is, where a is the 
Pauli spin operator and /is = 2^ • Here the spin of the 
electron can have two possible orientations one is j" or \. 
with respect to the applied magnetic field H. We assume 
the magnetic field is in the z direction, ie., H = H z . The 
Hamiltonian with the inclusion of spin can be obtained 
by adding — hb{&-H) to Eq.(49). The total Hamiltonian 
is given by 



2m 



(P 



-A) 5 
c 



1 22 

2 



(61) 



The first term in the Hamiltonian gives rise to diamag- 
netism and the last term is responsible for the paramag- 
netism. With Eq.(61), now the energy eigenvalues can be 



written as E nuri2yCr = E ni ^ 2 ■ 
partition function is given by 



/ib&zH, Hence the total 



F= I(ln(2sinh(^±))+ln(2sinh(^))) , (57) 



the internal energy is given by 



U 



flbj-l 



coth 



coth 



BHuj_ 



Entropy 
k B 



(58) 



ln(2sinh(^±))-ln(2sinh(^)) 



+ 



— - — coth( — - — ) + 



--— coth(--— ) (59) 



Z = Z + {fi,u + )Z_{p ) uJ)Z.{p) , (62) 
where Z S (B) is the factor introduced by the spin: 

Z S (B) = Tre'^-^ 11 ^ 

rp ( \ 

~ lr [ e~^ H ) 

= 2cosh(/3/i B #) . (63) 
Hence the total partition function turns out to be 



2cosh(^) 



2sinh(^±).2sinh(^; 



(64) 



The thermodynamic quantities: Free energy 
F 



I { ln(2sinh(^±)) 



ln(2sinh(^)) 



ln(2cosh(^))} 



(65) 



U 



- coth( 
tanh( 



2 

2 ' 



the internal energy 

tkj. 

~2 

hiu, 

Also entropy can be calculated as 

5 
k B 

f3hio + 
2 

2~ 

and the specific heat 



^coth( — 



(66) 



ln(2sinh(- — ±)) - ln(2sinh( P 



ln(2cosh(^))- 
— coth( — 



coth( 



2 
2 



)) 



tanh(^) ,(67) 



^ = (^) 2 cosech 2 ( — ) 

+ (^)Wh^) 

- (^) 2 sech 2 (^) . 



(68) 



Again in the limit of vanishing magnetic field, the re- 
suts for a free electron (in two-dimensions) in a magnetic 
field is recovered. 



IV. RESULTS AND DISCUSSION 
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FIG. 1: The variation of internal energy with respect to the 
temperature, for each oscillator. 
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FIG. 2: The variation of specific heat with respect to the 
temperature 



A. Bosonic, Fermionic and Supersymmetric 
Oscillators 

Bosonic oscillator: The behavior of thermodynamic 
properties of a simple bosonic harmonic oscillator both 
at low and high temperatures, is well studied. For a com- 
prehensive explanation we discuss these here. As T — > oo 
(ie., (3huj — > 0), the average energy U ps /c b T + clas- 
sical value is exactly recovered. In this problem we have 
a characteristic parameter , which is the ratio of the 
energy levels of the oscillator to the mean thermal en- 
ergy available. If fc^T is very large, discrete structure of 
the energy spectrum does not show up, one would ob- 
tain the classical value. In this limit, the specific heat 
converges, as expected. On the other hand, in the limit 
of T — > 0, we may see the largest deviations from the 
classical case in the internal energy and it turns out to 
be U ps 4p. Apparently, the specific heat vanishes like 
x 2 e~ x (x — > oo). 

The vanishing of specific heat at low temperatures is a 
characteristic feature of all quantum systems with a dis- 
crete energy spectrum. One way of explaining this effect 



is: if the thermal energy fc^T is much smaller than the 
oscillator energy fku, then an oscillator cannot (only with 
extremely small probability)be excited to the higher en- 
ergy level. Now the system is absolutely incapable of ab- 
sorbing energy provided by the heat bath. But classical 
systems can absorb arbitrarily small energies fc^T. Vari- 
ation of entropy is as expected. At low temperatures, 
the entropy approaches zero and at high temperatures 
entropy is infinite. 

Fermionic Oscillator: At low temperatures, the inter- 
nal energy of the system is like U ps ^f^, ie., the system 
tries to remain in the ground state. At high temper- 
atures, the variation of internal energy with respect to 
temperature is U ps gjjHf > means the syetm tries to pop- 
ulate the two available energy states. The specific heat of 
the system is vanishingly small at low temperatures, and 

that vary in a form which is given by 9 s - ps (^^) 2 e~ k B T . 
The specific heat vanishes at high temperature as well in 
a way ^ ps ( 2 ^" t ) 2 . Somewhere around T = r-, it dis- 
plays a maximum. Writing A for the energy difference 
between the two allowed states of the system, we may 
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FIG. 3: The variation of entropy with temperature 

write the specific heat as 

^ = (r^;) 2 e^(l + e^)- 2 . (69) 
kb k B I 

A specific heat of this form is generally known as the 
Schottky specific heat; characterized by an anomalous 
peak, and it is observed in all systems with an excita- 
tion gap A. The variation of entropy is as expected, it 
is vanishingly small at low temperaures (ksT << hu); 
it rises to a maximum rapidly when fc^T is of the order 
of the energy difference between the two states and it 
approaches a limiting value ks In 2 at ksT » Hlo. 

Supersymmetric oscillator: The variation of the inter- 
nal energy of a supersymmetric oscillator is almost like 
that of a bosonic oscillator. Here the specific heat is 
also vanishing exponentially at low temperatures just as 
in the case of a bosonic oscillator. Here also we can 
see at some temperature T, the specific heat displays 
a maximum. The specific heat of the supersymmetric 
oscillator vanishes exponentially at low temperatures as 
~ 2(-^Hp) 2 e k B T , ie., twice as that of a fermionic (or 
bosonic) oscillator at low-temperatures. At high temper- 
atures, the variation of specific heat is like (1 + jr^) 2 , 
which is the square of the respective quantity of a bosonic 
oscillator at high temperatures. The entropy vanishes ex- 
ponentially (« (3tkjje~P hu> ) at very low temperatures, and 
goes to infinity at very high temperatures. 

The variation of internal energy, entropy and specific 
heat of all the three oscillators with respect to tempera- 
ture are plotted in Fig.l, Fig. 2 and Fig. 3. 



B. Harmonically bound electron in a magnetic field 

We have two temperature scales, one is ksT/huj , 
which appears in the oscillator problem and the other 
is ksT/huJc, set by the magnetic field. It is to be noted 
here that the arguments of the hyperbolic terms in the 
expression for the specific heat contains both the temper- 
ature scales. We plot the thermodynamic quantities with 
respect to the temperature scale ksT '/Ituq. Let us define 



FIG. 4: The variation of specific heat when uu c = OAuq, with 
temperature 




FIG. 5: The variation of specific heat when uu c — 4u;q, with 
temperature 



— = and T = 
and we may write 



fceT 'jhujQ. So that fvjj c /2kBT = 



1 

wT 



If 



2 2T 

the cyclotron frequency is much greater than the os- 
cillator frequency ie., lo c » ujq, then oj + w uj c and 
u;_ sa 0. That means a magnetic dominancy is clear. 
At very low-temperatures, the internal energy goes as 

U ~ %\J cjq + ^p, and if the magnetic field is high, ie., 

lu c » u>o, then it is approximately equal to which 
is similiar to the one dimensional harmonic oscillator low- 
temperature expression with the oscillator frequency be- 
ing replaced by the cyclotron frequency. If the magnetic 
field is very low, ie., lo c « loq, then U ~ hu>o, the ground 
state energy of a 2D isotropic oscillator. In addition to 
this results if we take the electron spin into account, then 
in high magnetic fields, internal energy seems to be van- 
ishing!, and produces the value Hu)q at very low mag- 
netic fields. At high temperatures, the internal energy 
(both in the absence and presence of spin) varies like 
U ~ 2/cbT + .... which is the expected classical value. At 
high temperatures, the specific heat varies with temper- 
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ature like 

(t )spin=Q = 2 — (-; — 



.(2a 



O(^), (70) 



and similiarly 
(— 1 



k R l 



2o; c 2 ) 



G 



0( 



y2 - 



(71) 



Clearly, it is understood that the classical equipartition 
theorem is satisfied at very high temperatures, with a 
degree of freedom of 2 (since this is a two dimensional 
problem) . 

At low temperatures the behavior of the specific heat 
for a spinless electron in the presence of a parabolic po- 
tential and magnetic field is just identical to that of two 
independant isotropic simple harmonic oscillators as it is 
evident from the Eq.(60). As T — > 0, the specific heat 
vanishes like 



[G8feo+)(e- 



)] 2 + [(^_)(e- 



)] 2 . (72) 



While the heat capacity at low-temperature in the case 
where the electron possess the spin angular momentum 
of J has an additional term due to the spin contribu- 
tion. This also goes like w [(/3ftw c )(e _ ^ £ )] 2 . Hence the 
specific heat (with spin), decays as 



C v 



■^ ± )] 2 + [( / 3^_)(e- 



(73) 



In the limit ujq » uj c (ie., very low magnetic field 
strength) both the hybrid frequencies oj + and w_ ap- 
proximately equal to uiq, the oscillator characteristic fre- 
quency. Apparently, for very low magnetic fields (5^- — > 
0), the heat capacity for both cases follow the same form 
upto the leading order, which is given by 



Cv 
k B 



2( 



Hloq 
2k B T 



) 2 cosech 2 ( 



2k B T 



(74) 



which is nothing but the specific heat of a simple har- 
monic oscillator with the factor of 2 due to the extra de- 
gree of freedom. As T — > 0, the above expression vanishes 
exponentially as 2{j3hiOQe z^ 11 ) 2 . If the spin is present, 
then the specific heat goes (at low-temperatures) like 



C u 



) 2 -(Phw c e~^f 



(75) 



which shows a clear indication that the harmonic oscil- 
lator contribution dominates over the magnetic contri- 
bution. This result is expected too. When the mag- 
netic field strength is comparatively small, the contribu- 
tion from the spin angular momentum of the electron 
{/J'BO'zH) is nullified by the other part of the energy 



(which is harmonic in nature) ie., E„ 



Heat capac- 



ity versus temperature is plotted for both low and high 




FIG. 6: The variation of entropy when lo c = 4u;o, with tem- 
perature 



values of magnetic fields, in Fig. 4 and Fig. 5 respectively. 
From the Fig. 4 it is clear that at very low temperatures 
and low magnetic field strength (ie., u> c << ujq), the 
specific heat shows the behavior of an Einstein oscilla- 
tor, ie., an exponential suppression. Conversely, as the 
strength of the magnetic field is increased, the amalga- 
mated curves in Fig. 4 started bifurcating after a certain 
value of the temperature. For very high values of the 
magnetic field, ie., when uj c » loq, then lu+ rs lu c and 
u;_ w 0. The variation of specific heat at low tempera- 
tures and in high magnetic field resembles the varaiation 
of the specific heat (with temperature) of a harmonic 
oscillator at low-temperatures, but with the harmonic 
oscillator frequency being replaced by the cyclotron fre- 
quency, ie., 



k B 



(/3hui c 



(76) 



From Eq.(76), it is evident that ^ k at high mag- 
netic fields and at low-temperatures, and we can observe 
a small and sudden drop or dip in specific heat as in- 
dicated in Fig. 5. But neverthless, when the system is 
attenuated with the thermal fluctuations, ie., at high 
temperatures, all contributions are nullified and we can 
observe a clear resumption of the classical equipartiton 
value. The entropy variation is as expected. At very 
low-temperatures, the degree of disorder in the system 
is negligible and the system is said to be in the ground 
state. At low-temperatures the contributions from each 
terms of the Eq. (59) and (67) cancels each other to leave 
no contribution to the entropy by validating the third 
law of thermodynamics. We observe from the Figs. 6 and 
7, that both the curves follow same path at very low- 
temperatures and get bifurcated and approaches infinity 
at high temperatures. Variation of entropy with temper- 
ature is given in Fig. 6 and Fig. 7, both for high and low 
magnetic field strengths. 
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FIG. 7: The variation of entropy when u) c = (0.4)uo, with 
temperature 



heat and the entropy of these systems respectively with 
respect to temperature. All the low-temperature prop- 
erties are in conformity with the third law of thermody- 
namics and particulary specific heat showed an exponc- 
tial suppression in the low-temperatures. We have done 
a detailed calculation of the thermodynamic quantities of 
the well studied Landau-Fock-Darwin model which can 
be effectively converted into a harmonic oscillator prob- 
lem, both in the absence and presence of the spin angular 
momentum of the electron and the graphical results are 
in much coordination with the analytical results. We 
have discussed in detail both the low-temperature (both 
in the low (u> c << ujq) and high (u> c >> ujq) magnetic 
fields) and high temperature variations of the respective 
thermodynamic quantities. In view of the nanoscopic or 
mesoscopic systems, the LFD model has got numerous 
applications and is worth to study its thermodynamics. 



V. CONCLUSION 

We have examined the thermodynamics of the toy 
models of quantum mechanics viz., a bosonic harmonic 
oscillator, a fermionic oscillator and the combination of 
the both, a supersymmetric oscillator. We have graphi- 
cally analysed the variation of the internal energy, specific 
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